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ABSTRACT 

A  simple  test  is  given  for  determining  whether  a  given  matrix 
is  the  eigenvector  matrix  of  an  (unknown)  unreduced  symmetric 
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L  INTRODUCTION. 

In  a  sequence  of  papers  during  the  1960's  Slepian,  Landau,  and  Pollack 
investigated  functions  that  are  nearly  time  limited  ( supp  f  e[-7\7’])  and  nearly 
band  limited  ( supp  f  e(— 0,(1]),  see  [S.,1983].  Among  other  things  they 
discovered  that  the  integral  operator  K  on  L 8[— 7\7’]  given  by 

(A/XO  :=  flT  ~^7I/(s)  * 

commutes  with  a  second  order  differential  operator  D  given  by 

(!>/)<*)  :=  ((  7*-**)/ ’(*)]’  -  Cftr8/(*) 

and  zero  boundary  conditions.  Here  / 1  denotes  the  derivative  of  / .  Moreover  D 
has  a  simple  (point)  spectrum.  Consequently  the  eigenfunctions  of  K  are  the 
eigenfunctions  of  D  .  The  authors  managed  to  generalize  this  result  from  x  and 

t  in  R  to  *  and  t  in  RN.  A  decade  later  Slepian  demonstrated  the  discrete 
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analogue  [S.,1978]:  K  becomes  a  symmetric  Toeplitz  matrix  with  (i,j)  element 

,  -n  ^  n 

and  D  becomes  a  symmetric  tridiagonal  matrix.  This  is  a  rare  example  of  a  Toe¬ 
plitz  matrix  which  can  be  expressed  as  a  polynomial  in  a  non-Toeplitz  matrix. 

Recently  Grunbaum  has  been  extending  these  results  in  several  directions. 
For  example,  in  [G.,  1981]  he  describes  all  symmetric  Toeplitz  matrices  which 
admit  some  unreduced  tridiagonal  matrix  in  their  commutators.  It  turns  out 
that  each  such  matrix  is  determined  (essentially)  by  just  three  of  its  elements: 
(1,2).  (1,3),  and  (1.4).  In  [G.,  1982]  he  exhibits  Hankel  matrices  (including  the 
Hilbert  matrix)  that  commute  with  unreduced  tridiagonals. 

It  is  well  known  [M.  8c  M.,  p.77]  that  if  A  commutes  with  M  and  A  has  simple 
spectrum,  then  M  is  a  polynmial  in  A  and  they  share  the  same  eigenvectors. 
What  permits  a  symmetric  matrix  A  to  have  a  tridiagonal  in  its  commutator  is 
that  its  set  of  eigenvectors  is  very  special.  The  matrix  of  eigenvectors  may  be 
taken  as  orthogonal  and  this  communication  characterizes  the  eigenvector 
matrix  of  a  symmetric,  unreduced  tridiagonal  matrix  by  means  of  a  simple  test. 
It  turns  out  that  the  first  two  rows  (or  the  last  two  rows)  determine  all  the  oth¬ 
ers,  because  they  fix  the  eigenvalues. 

After  introducing  notation  and  reducing  the  general  problem,  we  give  a  list 
of  known  properties  of  these  eigenvectors.  Our  characterization  is  given  in  Sec¬ 
tion  5.  A  bound  on  the  number  of  zero  elements  in  a  row  is  given  in  Section  8 
along  with  a  result  connecting  eigenvalue  separation  and  the  top  row  of  the  nor¬ 
malized  eigenvector  matrix. 

&  NOTATION  AND  DEFINITIONS 

Capital  roman  letters  denote  matrices,  with  A  and  T  reserved  for  real  sym¬ 
metric  matrices.  The  elements  of  B  are  written  B(ij').  Lower  case  letters  x,  y. 


sin  2no(i-j 
2fr  (i-;) 


...  denote  column  vectors,  lower  case  Greek  letters  a,/S . denote  scalars,  but 

i.  j ,  k ,  l,  m,  n  are  reserved  for  indices.  G*  denotes  the  transpose  of  G. 


A  matrix  C  is  tri diagonal  if  C(i,j )  =  0  whenever  |i— j  |  >  1.  A  tridiagonal 
matrix  C  is  unreduced  if  C(i,  i+l)  *  0  and  C(i  +  1,  i)  *  0  for  all  relevant  indices 
i.  The  spectrum  of  C  is  the  set  of  eigenvalues  of  C,  ignoring  multiplicity. 

For  brevity  let  UST(n)  be  the  set  of  nxn  real  unreduced,  symmetric,  tridi¬ 
agonal  matrices.  Let  UST+(n)  denote  the  subset  of  UST  in  which  the  (i.i  +  1)  ele¬ 
ments  are  positive. 

Eigenvectors  are  defined  only  to  within  a  scalar  factor.  However,  the  eigen- 
vecrtorsj  of  a  real  symmetric  matrix  may  be  chosen  to  be  orthonormal,  and  this 
is  a  natural  convention  that  we  too  will  follow.  When  A  has  simple  spectrum  we 
may  speak  of  the  eigenvector  matrix  G  and  write 


A  =  GAG* 

where  A  is  diagonal  and  G  is  orthogonal.  For  a  given  ordering  of  the  eigenvalues 
in  A  the  matrix  G  =  (gj,  ...,.-gn)  is  unique  (up  to  replacement  of  gt  by  -gi  ).  If 
T  =  Ti>n  eUST  we  define  its  submatrix  Tj>n  by 


The  characteristic  polynomial  of  T^n  is 


ai  Pi 

| Pi  <*i+i  Pi+i 

Pi*\ 


'  Pn-l 

Pn-l  an 


Xi.n(0  :*  det[fI-Tf.n] 

is  abbreviated  by  x  or  yn.  The  spectral  factorization  of  T  =  Tl  n  e  UST(n)  is 
written 


T  =  S0S* 

where  0  :=  diag^x . i>n)  with  eigenvalue  ordering 

*i<*2<  ••  •  <fln, 

and  S  :=  (sj . sn)  is  the  orthogonal  eigenvector  matrix.  Thus 


4- 


Tst  =  Si  i>t,  i  =  1,  ....  n.  ands/sf  =  1. 

a  PROBLEM  REDUCTION. 

The  identity  matrix  I  is  tridiagonal  and  consequently  any  orthogonal  matrix 
is  am  eigenvector  matrix  of  some  tridiagonal.  This  case  shows  the  importance  of 
insisting  on  a  simple  spectrum  for  T. 

Consider  next  tridiagonal  matrices  T  which  are  reduced 
(tO'.J+1)  =  0  for  some  and  yet  have  simple  spectrum.  Such  Ts  are  the  direct 
sum  of  two  or  more  unreduced  tridiagonals,  say  T  =  diagiT^  T2>  Ta),  with  TteUST. 
Its  eigenvector  matrix  is  the  corresponding  direct  sum  of  the  eigenvector 
matrices  of  Tjt  Ta.  and  T^  providpd  that  the  eigenvalues  of  Tt  are  kept  together  in 
the  ordering  of  the  eigenvalues  of  T.  However,  no  matter  what  the  ordering  of 
the  eigenvalues,  Ts  eigenvector  matrix  is  not  fully  indecomposable.  We  follow 
the  terminology  of  [M.  8c  M.,  1964]  here;  some  authors  would  say  that  Ts  eigen¬ 
vector  matrix  is  reducible. 

The  procedure  for  deciding  whether  a  given  orthogonal  Q  is  the  eigenvector 
matrix  of  some  tridiagonal  with  simple  spectrum  is  to  column-permute  Q  into  a 
direct  sum  of  its  fully  indecomposable  parts  and  then  to  determine  whether 
each  part  is  the  eigenvector  matrix  of  some  Tt  in  UST. 

4.  KNOWN  PROPERTIES. 

FACT  1.  If  TeUST(n)  then  T  has  n  distinct  real  eigenvalues. 

Proof.  The  minor  of  the  (l,n)  element  of  T-fl  is  *fjf  T(t+,i)  which  vanishes 

<■1 

for  no  value  of  (.  So  rank  (T-fl)  k  n-1.  The  eigenspace  of  each  eigenvalue  has 
dimension  1.  The  result  follows  from  the  symmetry  of  T.  *. 

YARNING:  The  eigenvalues  of  T  may  be  very  close  indeed,  even  when  the  off 
diagonal  elements  are  not  small  relative  to  the  eigenvalue  spread,  but  see 


Theorem  3  in  Section  6. 


FACT  2.  If  TeUST(n)  then  the  eigenvalues  of  Ti>n_!  separate  those  of  Ti.n.  The 
eigenvalues  of  T2>n  separate  those  of  Ti,„.  ». 

Proof.  See  [  P.,  1980,  Chap.  10]. 

The  next  result,  which  uses  the  notation  established  in  Section  2,  was 
extended  by  Gantmacher  and  Krein  [G.  k  K.  1950]  to  the  class  of  oscillation 
matrices.  Karlin  has  extended  their  results  to  include  differential  operators  [K., 
19]. 

Recall  that  the  eigenvalues  of  T1-n  are  ordered  by 

v 

tfi  <1*2  <  •  •  •  <  d„. 

We  sketch  a  proof  because  it  is  short. 

FACT  3.  If  TeUST^Ti)  then  the  number  of  sign  changes  between  consecutive 
components  of  ^  ’s  eigenvector  s4  is  n-i,  for  i  =  1 . n. 

Proof.  By  the  3  term  recurrence  governing  the  Xi.t  one  can  verify  that  s*  is 
proportional  to 

[l.Xi(*i)/0i.  . •• 

and  each  ft  >  0.  The  polynomials  (l>Xi’--->Xn-ii  form  a  Strum  sequence.  (See 

[W..  1965,  p.  300]  but  observe  that  det[’T-('I]  =  (-l)nXn(f)  )•  It  follows  that  the 
number  of  sign  differences  in  1  l.Xi(0.-.Xn-i((")i  equals  the  number  of  eigen¬ 
values  of  Tl-n  greater  than  f  ,  not  less  than  f  .  The  result  follows  from  Fact  2. 

The  following  more  detailed  results  can  be  found  in  [P.t  1980,  Chap.  7]. 

FACT  4.  (C.C.  Paige).  If  TeUST(n)  then 

а)  *<(1)V(*<)  s  Xs.»(**). 

б)  *i(n)V(*<)  =Xn-i(^<). 

c)  *i(l)st(n)x’(^t)  =  Pi  fa  •  •  fa-i 

COROLLARY.  *(1)  #  0 ,St(n)  *  0„  i  =  1 . n. 

Proof.  By  Fact  2,  X23OM  *  0.  Xn-i(^t)  *  0.  By  Fact  1,  x'(^i)  *  0.  •. 


The  corollary  is  a  well  known  necessary  condition  for  an  orthogonal  matrix 
to  be  an  eigenvector  matrix  for  some  TeUST  .  In  the  next  section  we  generalize 
this  result  to  the  other  rows. 

If  TeUST(n}  then  there  is  no  loss  of  generality  in  assuming  that  TeUSTv(n). 
More  precisely,  if  TeUST  then  there  is  a  diagonal  matrix  A  with  A(i  ,i )  =  ±1  such 
that  ATA  €  UST*  .  If  S  is  the  eigenvector  matrix  of  T,  then  AS  is  the  eigenvector 

matrix  of  ATA.  Moreover,  we  may  take  1)  >  0.  i  =  1 . n. 

FACT  5.  (Uniqueness  of  Reduction.)  If  Q‘ AQ  =  T  e  UST+(n)  ,  and  Q  is  orthogonal, 
then  Q  and  T  are  determined  by  A  and  q  t  or  by  A  and  gn . 

Proof.  See  [P.,  1980,  Chap.  7], 

COROLLARY.  Let  S0S*  =  TeUST  (n)  be  the  spectral  factorization  of  T.  Both 
T  and  S  are  determined  by  0  and  e\ S  or  by  0  and  e£S.  Our  last  fact  makes  expli¬ 
cit  the  dependence  of  Q  on  gt  and  A  in  Fact  5.  The  formula  goes  back  at  least  to 
C.  Lanczos. 

FACT  6.  Let  T  =  Q*AQ  e  UST+(«)  with  Q  orthogonal. 

?j+i  =X/(A)gi/(ft  '  '  '  /*/)•  3=1 . n-1. 

Proof.  See  [P.,  1980,  p.  116]. 

5.  CHARACTERIZATION  OF  THE  EIGENVECTOR  MATRIX 

There  is  a  whole  family  of  tridiagonals  which  share  a  common  eigenvector 
matrix. 

LEMMA  1.  If  TeUST*  then  AT-aleUST*  for  all  a  and  all  A  >  0.  Moreover  T  and 
AT— crl  have  the  same  eigenvector  matrix. 

Proof.  The  first  asertion  is  trivial.  The  second  uses  the  orthogonality  of  the 
eigenvector  matrix  S;  T  =  S0S1  implies  AT-al  =  S(A0— oI)S* .  ■. 

This  trivial  lemma  is  the  key  to  a  simple  characterization  of  S.  It  shows  that 
in  the  search  for  a  T  there  is  no  loss  in  taking  aj  =  T(l,l)  =  0  and  /?j  =  T(l,2)  =  1. 


THEOREM  1.  Let  G  =  (gi . ....  gn)  be  orthogonaL  G*  is  the  eigenvector  matrix  of 
some  TeUST(n)  if  and  only  if  the  following  conditions  hold. 

(1)  gi(i)  *  0.  i  =  1 . n. 

(2)  The  numbers  iS*  are  all  distinct,  where  :=  g  2(i)/  g  i(i).  i  =  1 . n. 

(3)  G*0GeUST,  where  9  =  diagi&i,  ....  $n). 

Proof.  Condition  1  was  established  as  the  Corollary  of  Fact  4. 

Because  G  is  orthogonal,  G4  9G  =  T  can  be  rewritten  0G  =  GT.  Equating 
column  1  on  each  side  yields 

901  =  01<*1  +0201. 

By  Lemma  1  there  is  no  loss  in  taking  at  =  0,  /9j  =  1,  to  obtain 

=  02^).  i  -  1,2, . n.  Thus  0  is  determined  by  gY  and  gz.  By  Fact  1  the 

must  be  distinct  This  is  condition  2. 

If  G4AGeUST(n)  for  some  A  =  diag(\u  ....  **),  then,  by  the  previous  para¬ 
graph,  A0i  -0 i«*i  +02^i  for  some  ai.jSj  and  thus  \  ~  at  +  ,i  =  1 . to. 

Consequently  G*0G  =  (G4  AG— a^)/  /?i  ellST(n )  as  well.  This  establishes  condition 

3.  * 

Condition  3  is  easily  stated  but  will  fail  for  most  orthogonal  matrices  satis¬ 
fying  conditions  1  and  2.  We  can  replace  (3)  with  a  sequence  of  necessary  condi- 
tons.  As  soon  as  one  fails,  we  know  that  G(  cannot  be  an  eigenvector  matrix. 

DECISION  PROCEDURE.  Given  distinct  values  i5t  :=  02(i)/0i(i),  i  =  1 . n, 

and  0  -  (Hag (tfj, ....  tfn)  proceed  as  follows: 

for  j  =  2 . n— 1,  do 

formjfj+i  :=  0p,  - g tat  -0*- iff/.*, 
where  a t  :=  gj&gJw  /?/.  a/.  i}, 

If  9j*\  *  ±0/+ \/  Il0>+i  ||  then  exit  and  report  failure, 
if  no  failure  then  G*  is  an  eigenvector  matrix  for  the  T  defined  by  the  a*  and  /Jt. 


Proof.  By  Fact  5  (uniqueness  of  reduction)  there  is  a  unique  orthogonal 
matrix  G  and  a  unique  T+  in  UST+  satisfying 

T+  =  G‘0G, 

All  T  in  UST  orthogonally  similar  to  3  sure  given  by 

T  =  AT+  A  =  AG*  QGA 

with  diagonal  A  and  A(i,i)  =  ±1.  Hence  G*  QGeUST  if  and  only  if  G  =  GA  for  some 
A. 

The  algorithm  in  the  procedure  is  simply  the  construction  of  the  matrix  G 
column  by  column  together  with  the  test  gt  =  ±gt.  «. 

There  is  a  result  dual  to  Theorem  1  which  uses  the  last  row  instead  of  the 
first  one.  We  state  it  without  proof. 

Theorem  1’.  Let  G  =  (pi,  ....  gn)  be  orthogonal.  G*  is  the  eigenvector  matrix 
of  some  TeUST(n)  if  and  only  if  the  following  conditions  hold. 

V  9n(i)*  0.i  =  l . n. 

2. ’  The  numbers  ?*(:=  ffn-i(0/ffn(i))  must  be  distinct. 

3. ’  G*$GeUST,  where  $  =  diagfo^  ....  <pn). 

Since  different  normalizations  are  involved,  we  will  not  have  ^  =  <pi  in  gen¬ 
eral.  However  they  must  be  lineraly  related, 

^  -  fin- i9i  =  «Xn,  i  =  1,  ....  n. 

In  practice  the  decision  procedure  could  work  simultaneously  from  the  top 
row  and  the  bottow  row  towards  the  middle. 

ft  ZERO  ELEMENTS. 

THEOREM  2.  If  S  is  the  eigenvector  matrix  of  TeUSF(n)  then  the  number  of 
zero  elements  in  its  j’th  row  is  at  most  min  (j-1,  n-j). 

Proof.  Since  T  ~  AT  A  for  T+cUST*  and  A(i,i)  =  ±1  there  is  no  loss  in  consid¬ 
ering  TeUST+(n).  Apply  Fact  8  to  the  relation  T  =  S9S*  to  conclude  that 


S*8;>iO?i  •  '  '  Pj)  =  Xj(Q)Stel,  i.e., 

s<0  +  l)(/?i  •  •  /?/)=  XjCtfiKC1).  *  -  1 . n. 

Now,  is  a  polynomial  of  degree  j  and  independent  of  tV  It  can  have  at  most  j 

zeros  and  so  S(J  +  1,.)  can  have  at  most  j  zero  elements.  These  occur  when,  and 

only  when,  an  eigenvalue  of  T  is  also  an  eigenvalue  of  a  submatrix  T]  j . 

By  Fact  5  (uniqueness  of  reduction)  T  and  S  are  also  determined  by  0  and 
sn .  The  analogous  formula  is 

Si(j  +  l)(pj+i  ■  ■  ■  Pn)  =  Xj+1.n(^i)Si(n).  i  =  1 . n. 

Since  X/+i,n  is  of  degree  n-j-1,  S (J  +  1,.)  can  have  at  most  n-j-1  zero  elements  and 

these  occur  when  an  eigenvalue  of  T  is  also  an  eigenvalue  of  T^+u».  •. 

Our  last  result  concerns  the  separation  of  the  tf’s . 

THEOREM  3.  For  T  eUST+(n), 

ft  ~*i)/  Pi  =  1/  fts(l,i  )  fe  n* 

t.i*l  »*! 

Proof.  Define 


B\-dxag{\,Pi,PiPz....,PiPz..Pn-i)> 

E:=cfiag'(S(l,l),S(l,2),  •  •  •  ,S(l,n)). 

The  displayed  equation  in  the  proof  of  Theorem  2  implies  that 

S=Z?-1CE,  where  C(i,jO  :=  Xi-iCfy)- 
By  standard  properties  of  determinants  and  Vandermondes 

detC  =  detV  =  ft  (*fy- iJt)  >  0, 

iJ-l 

<<S 

where  V(i  Since  S  is  orthogonal, 

±1  =  detS  =  detC  x  detS/detff. 

By  our  normalizations,  including  S(1,Z )  >  0,  the  right  side  is  positive  and  so  det  S 
=  +1.  Finally 


deti?  =  pr'Pr 2  ■■■Pn- 1  =  ft  /Si 


The  last  inequality  comes  from  the  generalized  Arithmetic-Geometric  mean  rela¬ 
tion.  • 

We  quote  without  proof  the  analogous  result  using  the  last  row  of  S. 
THEOREM  4.  For  T  €USJ\.(n), 

ft  =  l/ft|S(n.Z)|  few* 

i.i=l  1=1 

<<i 

These  results  are  informative.  As  so  can  one  or  more  of  the 

differences  ,di+1wdi-*0  without  violating  the  lower  bound.  Of  more  interest  is  the 
case  of  known  fixed  fts.  The  differences  of  tf’s  is  bounded  above  by  the  spread 
i)n— which  is  itself  bounded  by  2||T|].  This  yields  a  crude  but  explicit  lower 
bound  on  minf'ft^i^ft)  in  terms  of  the  fts  and  ||  T  ||.  Although  eigenvalues  of  T*s 
in  UST(7i)  can  be  surprisingly  close  to  each  other,  they  cannot  be  arbitrarily 
close  when  the  ft s  are  given. 

The  relationships  are  even  clearer  in  the  special  but  important  case  of 
ft  =  1  for  all  i.  For  small  values  of  n  the  presence  of  tiny  S(l.f)  precludes  close 
t?’s  but  as  n  increases  this  effect  weakens  rapidly.  Wilkinson’s  matrix  W21  is  an 
example  of  this  phenomenon.  See  [W,p.308]  for  more  information  on  this 
matrix. 

References 

D.  Slepian  (1983), "Some  Comments  on  Fourier  Analysis,  Uncertainty  and 
Modeling", S.I.A.M.  Rev.  25  (1983)  379-394. 


D.  Slepian  (1978),  "Prolate  Sheroidal  Wave  Functions,  Fourier  Analysis,  and 
Uncertainty",  Part  V,  Ibid,  vol.  57,  pp.  1371-1430. 


F.A.  Grunbaum  (1981),  "Eigenvectors  of  a  Toeplitz  matrix:  discrete  version  of  the 
prolate  spheroidal  wave  functions",  SIAM  J.  Aig.  Disc.  Math  2,  pp.  138-141. 


F.A.  Grunbaum  (1981),  "Toeplitz  Matrices  Commuting  with  Tridiagonal  Matrices", 
Lin.  Alg.  Appl.  40,  pp.  25-36. 

F.A.  Grunbaum  (1982),  "A  Remark  on  Hilbert’s  Matrix",  Lin.  Alg.  Appl.  43,  pp. 
119-124. 

B.N.  Parlett  (1980),  The  Symmetric  Eigenvalue  Problem  ,  Prentice  Hall  Inc.,  New 
Jersey. 

J.H.  Wilkinson  (1965),  The  Algebraic  Eigenvalue  Problem.  Oxford  Univ.  Press. 

F.P.  Gantmacher  and  M.G.  Krein  (1950),  Oscillation  Matrices  and  Kernels  United 
States  Atomic  Energy  Commission,  Office  of  Technical  Information,  AEC- 
tr-4481  Physics  (1961).  The  Russian  original  appeared  in  1950. 

M.  Marcus  and  H.  Mine,  A  Survey  of  Matrix  Theory  and  Matrix  Inequalities.  Allyn 
and  Bacon,  (1964). 


*■ 


0 


